Strong constraint on dark matter — neutrino interactions via light vector bosons 
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Dark matter velocity-dependent self-interactions and interactions with neutrinos mediated by a 
new MeV-scale vector boson V have been shown to solve all the small-scale-structure problems 
in a ACDM cosmology [van den Aarssen et al. [1]; Phys. Rev. Lett. 109, 231301 (2012)]. We 
show that in the minimal version of the model such interactions with Standard Model neutrinos 
lead to 'V-strahlung" from neutrinos in the final states of kaon and W decays, producing vastly 
too many missing-energy events. The minimal formulation of the proposed model is thus strongly 
excluded. More generally, we present new constraints on the coupling strength to neutrinos, g,^ < 
8 X lO^'' (mv/MeV), for a light vector boson of mass niv- These model-independent constraints 
provide one of the most stringent bounds on additional neutrino interactions. 



Introduction. — Precision measurements of the 
cosmic microwave background provide overwhelming 
evidence for dark matter (DM) as the dominant form 
of matter in the Universe [2-4]. These and other 
measurements at large distance scales are in remarkable 
agreement with the predictions of the Lambda Cold Dark 
Matter (ACDM) model of the Universe [5-7]. 

However, at the scales of galaxy clusters, galaxies, and 
yet smaller objects, ACDM predictions do not match 
the observations [8]. There are three important and 
enduring problems at small scales. First, "core vs. cusp" 
- flat cores are observed in the density profiles of dwarf 
galaxies, whereas numerical simulations predict sharp 
cusps [9-16]. Second, "too big to fail" - the most massive 
subhalos found in numerical simulations are denser than 
the visible subhalos of the Milky Way [17, 18]. Third, 
"missing satellites" - fewer satellite galaxies are observed 
than predicted in numerical simulations [19-27]. It 
has proven difficult to provide a solution - whether 
by using baryonic physics [28-36] or new particle 
physics [37-58] - to all three of these small-scale problems 
simultaneously while remaining consistent with the large- 
scale observations of ACDM. 

Recently, van den Aarssen, Bringmann, and 
Pfrommer [1] proposed a phenomenological model 
of DM, described by £ D —9xX^^X ^ g^i^V'^, that, 
remarkably, solves all of these small-scale structure 
problems. The key ingredients of this model are DM (x) 
with velocity-dependent self-scattering mediated by a 
new MeV-scale vector boson {V) and DM scattering with 
neutrinos mediated by V. Velocity-dependent XX ~^ XX 
scattering reduces the central densities of galaxies and 
their subhalos, addressing the "core vs. cusp" and the 
"too big to fail" problems. The X^ ~^ X^ scattering 
ensures the late kinetic decoupling of DM, addressing 
the missing satellites problem by suppressing structure 
formation for masses smaller than ~ 10^ Mq. The last 
requirement restricts the model to the parameter space 
shown by the shaded green region [1] in Fig. 1. 



Given the importance of the tension between the 
ACDM model and observations on small scales, it is 
urgent to test this possible resolution. However, this is 
quite challenging because the only other particles whose 
phenomenology is affected are neutrinos (in the model 
of ReL [1], Standard Model neutrinos; an extension to 
sterile neutrinos [1] is discussed below). We focus our 
attention on the neutrino interactions with V, both 




FIG. 1. Constraints on the boson mass my and the neutrino- 
boson coupling gi, from kaon decay (dashed red) and W decay 
(solid blue), as derived in this work. The hashed regions are 
excluded at > 90% C.L. The shaded green region is the now- 
excluded parameter space of the minimal version of the model 
in Ref. [1] that solves the missing satellites problem of ACDM. 



for testing this model and for developing more general 
constraints. Our results are shown in Fig. 1. 

Extra neutrino interactions. — It is interesting to 
ask if neutrinos have interactions besides those due to 
the weak force. Neutrinos (or dark matter) may be the 
best particles to use when searching for new interactions, 
because the electromagnetic and strong forces are absent. 
It is hard to detect neutrinos, but sometimes sensitive 
tests are possible anyway. 

In typical models of new interactions, including 
Ref. [1], it is assumed that the new bosons couple only 
to neutrinos, and not to the charged-lepton partners in 
their weak doublets, thus breaking SU(2)lx U(1)y gauge 
invariance explicitly. While somewhat unnatural (see the 
discussion in Ref. [59]), this can be done, and is necessary 
to avoid strong constraints from precise measurements of 
charged leptons. 

If neutrinos couple to new heavy bosons, this more or 
less just modifies the Fermi constant. However, a rich 
phenomenology is possible if neutrinos couple to new 
light bosons that are kinematically accessible. Models 
with a new light scalar boson coupled to neutrinos, 
e.g., a Majoron, have been extensively studied, and 
there are strong constraints on such couplings [60-76]. 
Interestingly, models with a new light vector boson seem 
to have been largely overlooked (except for Ref. [64], 
discussed below), which we rectify. 

A light vector boson instead of a scalar is employed in 
Ref. [1] to enhance the neutrino cross section with dark 
matter; the non-zero mass prevents 1/r^-type forces [77]. 

V can typically acquire a mass via spontaneous symmetry 
breaking. If V is an abelian gauge boson, the Stuckelberg 
mechanism could also provide a mass, but in that case 

V would, in general, kinetically mix with photons and 
can be probed using a variety of experiments. Here, 
we concentrate on the more general scenario of a broken 
gauge theory. 

An obvious way to constrain this model is to search for 
neutrinos from dark matter annihilation, but current and 
projected sensitivities [78-82] are not strong enough [1]. 
Instead, we show that it is possible to radiate bosons 
from final-state neutrinos in weak decays. This can turn 
a two-body process with a mono-energetic charged lepton 
into a three-body process in which the charged lepton has 
a continuous spectrum, indicating the presence of a new 
non-interacting particle to carry the missing energy. In 
addition, the total decay width will be increased. 

For a scalar boson, the rate of boson emission scales 
as ^ g^ , where g^, is the neutrino-boson coupling. 
In contrast, for a vector boson, the rate scales as 
~ gl{Edcca.y/mv)'^, where i^dccay is the available energy 
and my is the boson mass. This enhancement factor, 
which arises from the longitudinal polarization modes of 
the light vector boson in a spontaneously broken gauge 
theory, greatly increases the sensitivity to the coupling 
gi, (this enhancement is only possible due to the explicit 
breaking of Standard Model gauge invariance by coupling 
the V only to the neutrino [59]). We take the neutrino-^ 



coupling to be of the V — A form; our results have equal 
contributions from the vector and axial couplings. 

Constraint from kaon decay. — Kaons dominantly 
decay (branching ratio ~ 65%) via the 2-body leptonic 
channel K~ — ^ ^~I7^, for which the muon energy 
spectrum is a delta function in the kaon rest frame. 
If a new vector boson couples to neutrinos as g^V^^v, 
then there can be y-boson bremsstrahlung from the 
antineutrino in the final state if mv ^ 388 MeV; the 3- 
body decay K~ — )■ fjT v^V , has a dramatically different 
muon spectrum [83]. Because the muon spectrum from 
kaon decay has been measured very precisely, we can use 
this to derive a strong upper limit on the probability for 
V emission and hence also on g^- 

The setup of our problem is such that it satisfies 
all five criteria for the application of the narrow-width 
approximation [84], so the V is treated as an on-shell 
particle (it subsequently decays to vv). We therefore 
consider the 3-body decay K~ — ^ //^ F^ V ^ as shown in 
Fig. 2. Much of the calculation is similar to that for 
a related limit on parity- violating muonic forces [62]. 
Motivated by the model of Ref. [1], we develop new limits 
on general neutrino interactions. 



K- 




FIG. 2. Feynman diagram for K~ {us) decay to a muon where 
a 1/ is radiated from the final state antineutrino. The hadronic 
matrix element {0Jm7"(1 — 75)s|_K'~) = fKPK is denoted by 
the shaded circle. 

The amplitude-squared for this process, after summing 
over final spins, is 



|X|2 



TO^ (Wy +2p„ -PV)^ 

^{Pv-Pv) C^Pk -Pf^PK -P^ -mjfPp.-p^) 
-Amyp^ -pv {2{pK ■Pf.f -mxiPK • Pp +™m)} 
-my{2pK ■P^.Pk ■Pu-m^KPt.-Pi') , (1) 



where Pi and rrii denote the 4-momentum and mass of a 
particle i, Gp is the Fermi coupling, 0c is the Cabibbo 
angle, and /k is related to the hadronic matrix element 
given in the caption of Fig. 2. 

The terms are organized by powers of my- This 
is useful as my is smaller than all other masses in 
the problem. The largest contribution comes from the 



term that is proportional to l/rriy, as we have verified 
from the fuU calculation. This term arises due to 
the sum over polarizations of the V^-boson, which is 
~9fj.i^ + k^k^/rriy [59]. This polarization sum arises for 
a spontaneously broken gauge theory. The dependence 
of the decay rate on l/rriy is similar to that of the 
weak boson emission process t -^ bW, which scales as 
l/m2 



w 



The double-differential decay rate [86] in terms of the 
{Pi^ +PvY and m^g = (p^ +PvY 



Dalitz variables m 
is 



dT{K- 
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We integrate this over the allowed range of m.231 given 
in Eq. (40.22) of Ref. [86], to obtain the single-differential 
decay rate in terms of 771^2- Using the relation 771^2 = 
m\ -\- mi — 2mKE^, in terms of the muon total energy 
E^ in the kaon rest frame, the differential rate is 
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where N = g^ Gj, sin^ Be Ik /{^^Gt:-^ mj^). The main 
scaling of this result is ~ gt/fny, which means that our 
constraint in the {g^,mv) plane will be linear. 

In Fig. 3, we show the muon spectra from kaon decay in 
two cases: when V emission is forbidden (K~ — ^ /7~^m) 
and when it is allowed {K~ — > /i^T/^F). In both cases, we 
plot dV/dE^ normalized by the total (all modes) decay 
width Ftot- For the 2-body decay, the muons have a 
monoenergetic spectrum with E^ — 258 MeV; we show 
the measured result (including energy resolution) [83] 
in solid blue. For the 3-body decay, the muons have 
a continuum spectrum; we show this in dashed red, for 
f/i/ = 10~^ and 77iy = 0.5 MeV, which are a representative 
set of values favored in Ref. [1]. This produces events 
at energies where no excess events above the Standard 
Model background were observed (shaded region) [87]. 
We also show the approximate upper limit that we 
derived using a flat spectrum (in the energy range used 
for the search) from the upper limit as presented in 
Ref. [87]. 

To obtain our constraint, we use the results from 
a search for missing-energy events that was reported 
for kaon decays with muons observed with kinetic 
energies from 60 MeV to 100 MeV (£'^ between 
165.5 MeV and 205.5 MeV). We integrate our calculated 
differential decay rate, dT/dE^, over this range of i?^ to 
obtain the partial decay width F(A'^ — > pi^Vf^V). 
The constraint on the branching ratio 
r{K-^ IJ,- invisible) /F (iiT- ^ ^- z7^) < 3.5 x 10"*^ [87] 
leads to the limit on g^ shown in Fig. 1 (dashed red 
line). The constraint is 5^ < 8 x 10~^ (777^ /MeV) for 
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FIG. 3. Muon spectra from kaon decay for the standard 2- 
body decay K~ — >■ fi~ I?^ (solid blue) measured in [83] and the 
hypothetical 3-body decay K~ — > /7^ 17^ V (dashed red) with 
g^ = 10~^ and mv = 0.5 MeV. The shaded region shows the 
search region of Ref. [87], where no excess events were found. 
From this we derive an upper bound on the 3-body differential 
decay rate that is ~10^ times lower than the dashed red line. 



my ^ 388 MeV. This simple yet stringent constraint 
from kaon decay rules out the entire parameter space of 
Ref. [1] (shaded green region in Fig. 1). Only very low 
777y values, beyond the range shown in Ref. [1] and in 
our Fig. 1, may survive. 

This argument nominally only constrains the coupling 
of V to Vf^. However, if the neutrino- V^ coupling is 
independent of the neutrino flavor, as in Ref. [1], the 
limit applies to all three neutrino flavors. In fact, unless 
there is a special connection between flavor symmetry 
and dark matter, this should be true for any model. 

Constraint from W decay. — Our constraint can 
be made more robust by also considering W decay, 
which applies for all three flavor of neutrino and much 
larger values of 777y. In the following, we focus mostly 
on the places where there are differences from the 
kaon calculation. Similar considerations have been 
applied for electroweak bremsstrahlung in dark matter 
annihilation [88-90]. Our limits on general neutrino 
interactions are new. 

In the leptonic decay of the W boson, W^ — >■ £" Vi 
(branching ratio averaged over all three flavors ~ 10%) , a 
V boson can be emitted from the final state antineutrino, 
if a g^v^v coupling is allowed. The 3-body decay 
of the W-hoson (shown in Fig. 4) leads to additional 
events with missing energy, increasing the total decay 
width of the W . While this has been measured with 



only ^1% precision, the 771 
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leads to a comparable constraint on g^ for an MeV- 
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FIG. 4. Feynman diagram for W decay to leptons in which a 
V is radiated from the final state antineutrino. 



mass V (it may be possible to develop even more 
stringent constraints using the charged lepton spectrum). 
Although mw/mv ^ 1, the decay rate can be calculated 
perturbatively for g,^ <^ 1. 

The amplitude-squared for this process, after averaging 
over the initial spin and summing over final spins, is 
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The double-differential decay width [86] in terms of the 
Dalitz variables m^2 = {Pv + PvY and tti^s = (pi +pvY 
is 



dT{W~ -^ r Vn V) = 
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We integrate this over the allowed range of m\^ and also 
over the allowed range of to^2 to obtain the partial decay 
width V{W — ?> (.^viV). For the electron channel, this 
can be written as 
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to leading order in mv/'niwr. The width of the 2- 
body mode is T{W~ — >■ e~Pe) — G_f m^/(6'\/27r), so 
the relative rate of the three-body mode to the two- 
body mode is ~ g1{Mw /mvY / {Mtt'^). The factor 
(Mvy/my)^ allows us to probe very small values of g^. 

The experimentally-measured total decay width of the 
W is 2.085 ± 0.042 GeV [86], which agrees very well with 
the theoretically-calculated value, 2.091±0.002GeV [86]. 
If the rate of "V-strahlung" were too large, then 
the increase in the calculated total width would be 
inconsistent with experiment. To obtain a one-sided 90% 
C.L. upper limit on the neutrino-boson coupling y,^, we 
demand that r(T^- -^ f Vi V) < 1.28 x 0.042 GeV. The 



constraint is 5;^ ^ 17 x 10~^ (my/MeV) for mv ^ Ww 
and is shown by the solid blue line in Fig. 1. The 
constraint on g^ is approximately the same for all three 
lepton flavors. As above, this rules out the parameter 
space of Ref. [1], and only very low values of my may 
survive. A related constraint from Z decay is weaker due 
to destructive interference between the two contributing 
Feynman diagrams. 

Our constraint was derived assuming that single- 
V emission could be treated perturbatively. At the 
boundary we define, this is reasonable because the ratio 
of the leptonic decay width of the 3-body mode to the 
2-body mode is ~ 20% (and the ratio to the total decay 
width of the W^ is ~ 2%) and non-perturbative effects do 
not set in. To the right of our constraint in Fig. 1, this 
approximation will no longer be valid and the cascade 
emission of multiple V bosons will occur, for which non- 
perturbative methods must be used [59, 91]. In this 
region, which is still excluded, the decay rate will just 
grow more slowly than ~ g'^m'^/my. 

Conclusions and implications. — We present new 
and simple constraints on neutrino interactions with a 
light vector boson from the nonobservation of excess 
missing-energy events in the decays of kaons and W 
bosons. To the best of our knowledge, these are by far 
the most stringent constraints of this type. These model- 
independent constraints have a strong impact on the 
viability of models that make use of additional neutrino 
interactions. 

An interesting class of such models, proposed recently 
in Ref. [1], uses a vector boson to couple neutrinos to 
DM, in order to delay the DM kinetic decoupling and to 
provide a natural and elegant particle physics solution to 
the missing-satellites problem of ACDM. We rule out all 
of the interesting parameter space of the minimal version 
of the DM model in Ref. [1] . This opens up again the 
important problem of providing a unified solution to the 
small-scale problems of the ACDM model. 

Our constraints provide valuable guidance for future 
model building. If longitudinal modes of V are absent, 
these constraints are lifted. Our constraints can also 
be evaded if the V boson couples primarily to sterile 
neutrinos. However, the scattering of dark matter with 
neutrinos assumed in Ref. [1] would then require a large 
background of sterile neutrinos (or a larger g^,) and 
must respect strong limits on extra light particles [2]. 
If active-to-sterile transformations are efficient in the 
early universe, they might also be in supernovae, so 
limits based on neutrino scattering en route could be 
relevant [64]. The limit from SN 1987A is weaker than 
ours by several orders of magnitude, but the sensitivity 
could be greatly increased once supernova at cosmic 
distances [92] are detected. 
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